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S We present a simpler proof of the existence of an exact number of one 
i or more limit cycles to the Lienard system x = y — F (x), y = — g (x), under 

weaker conditions on the odd functions F (x) and g (x) as compared to those 
t-H available in literature. We also give improved estimates of amplitudes of the 

^ limit cycle of the Van Der Pol equation for various values of the nonlinearity 

s ! parameter. Moreover, the amplitude is shown to be independent of the 

asymptotic nature of F as \x\ — > oo. 
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1 Introduction 

There has been a considerable interest in the study of the number and nature of 
limit cycles in a Lienard equation 

x + f (x) x + g (x) = (1) 

recently [[7] — [IB]]- Limit cycles are isolated periodic curves in the phase plane 
and arise in numerous applications as self-sustained oscillations which exist even 
in the absence of external periodic forcing. The equation (TlJ is usually studied as 
an autonomous system, called the Lienard system, given by 

x = V ~ F ( x ) , V = -9 ( x ) ( 2 ) 
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where F (x) = j f (u) du. The phase plane defined by ^ is called the Lienard 
o 

plane. Lienard gave a criterion for the uniqueness of periodic cycles for a general 
class of equations when F (x) is an odd function and satisfies a monotonicity con- 
dition as x — > oo. An interesting problem for the system ^ is the determination 
of the number of limit cycles for a given odd degree (m) polynomial F (x). Lins, 
Pugh and de Melo [6J conjectured that the system ^ has at most N limit cycles 
if m = 2N + 1 or m = 2N + 2. Currently this problem is being investigated by 
many authors in connection with the still unsolved Hilbert's 16th problem. 

Giacomini and Neukirch |10) have developed a general procedure for construct- 
ing a sequence of polynomials whose roots of odd multiplicity are related to the 
number and location of the limit cycles of equation ([I]) when / (x) is an even degree 
polynomial. They have also given a sequence of algebraic approximations to the 
equation of each such cycles, although their method is mainly of experimental (nu- 
merical) in nature and a rigorous justification is still lacking. Hoist and Sundberg 
[T5] have extended Rychkov's theorem [5J for a class of F (x) having 5th degree 
polynomial like behaviour. The proof of Richkov's theorem however requires bi- 
furcation theory. Odani gave a proof on the existence of exactly N limit cycles 
of the Lienard equation ([!]) with g (x) = x. His proof does not make use of the 
bifurcation theory. His method also gave an improved estimate of the amplitude of 
a limit cycle. Recently there have been some progress in elucidating sufficient con- 
ditions extending the previous results. Chen and Chen [12] , for instance, proved 
the Lins-Pugh-de Melo conjecture for Lienard system with function F odd. On 
the other hand, it has been shown [16] that for suitable polynomial F of degree 
7, the system ^ has 4 limit cycles, contradicting the conjecture in [BJ. Chen, 
Llibrc and Zhang Q~7] proved a sufficient condition for existence of exactly N limit 
cycles for the system ^ with a general class of F (x) functions. We investigate an 
equivalent problem covering, however, a different class of functions F as compared 
to pi]. 

In the study on the number of limit cycles several authors have studied the 
equation ([I]) in the usual phase plane (viz. Theorem 7.10 — 7.12, Chapter 4 in 0]) 
whereas some considered the Lienard plane [IT]. In Theorem 7.10, Chapter 4 [3] 
the function / is taken as a periodic function. Theorem 7.11 is a generalization 
of Theorem 7.10 in which the function F' (x) = / (x) is a monotone function in 
certain regions. However Theorem [3] and Theorem [4] in the present paper do not 
depend upon the monotonicity of /. Rather, we have used the monotonicity of 
F. As a consequence, merely the sign of the function / determines the monotonic 
nature of F, and hence determines the number of limit cycles in Lienard system 
Thus our results cover a different class of functions than those covered by 
Theorems 7.10 and 7.11 mentioned above. Theorem 7.12, Chapter 4 in [4] and 
the theorem in [T7] have been proved on Lienard plane. Both of these results have 
assumed the existence of /3j £ [aj, aj+i], j = 2, 3, 4, . . . such that F (/3y) = F (ij-i) 
where, Oj's are positive roots of F and Lj's are unique extremum of F in [ay, aj+i] 
for j = 1,2,3, . . .. However, if we do not get any such (3j then these results are 
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not applicable. In such situations Theorem [4] in Section [4] of the present paper is 
still applicable to determine the exact number of limit cycles. One such example 
is given in section [5] 

In this paper we first give a simple but, nevertheless, an important extension of 
the Lienard's theorem for the unique limit cycle by removing the unbounded nature 
of the function F as x — > oo. Next, in Theorem [3] we prove that the system ([2| has 
exactly two limit cycles when the odd function F (x) undergoes two sign changes 
in x > and is monotonic not only as x —> oo, but also near (actually at the right 
of) the first zero. However, g (x) (g (x) > for x > 0) can be any odd continuous 
function. Example [6] in support of Theorem [3] reveals clearly the strength of this 
theorem over analogous results (e.g. Theorem 5.1) of pQ. The new insights gained 
from Theorem [3] (and also from Theorem [2| then provide a general approach in 
obtaining an existence theorem for multiple limit cycles in a systematic manner. 
In Theorem [4] we state a set of such conditions for the existence of exactly N limit 
cycles. Although we are dealing with odd functions F only, there are certain odd 
functions as shown in Example [7J which satisfy the conditions of Theorem [4] in the 
current paper but do not satisfy the theorem of [17] . This establishes our claim 
that the present theorems cover different classes of functions F than those covered 
in [4] and [17]. Moreover, as stated above, g (x) here is an odd function while for 
the theorem of [TTj g (x) — x. The second important result that we find in section 
an efficient upper estimate of the amplitude of the limit cycle for the system 

. The values of the amplitudes for the Van der Pol equation are obtained in 
Example [lj which are much more accurate compared to those in |llj and |14| . 

The paper is organized as follows. In section[2j we sketch the main steps of the 
proof of the classical Lienard theorem thus introducing our notations. In section 
[3j we discuss some special observations leading to an extension of the classical 
Lienard theorem. Our main result, Theorem [3j on the existence of two limit 
cycles is proved in section [4] In Theorem [4] we state the sufficient conditions for 
existence of exactly N limit cycles. An outline of the proof is given in Appendix 
(For a detailed proof, see [IS]). The proof of this general existence theorem is 
based on an induction method with non-trivial initial hypotheses for N = 1 and 
N = 2. We present some examples in section [5] highlighting the key features of 
the above theorems. Section [6] contains some concluding remarks. 

2 Lienard's Theorem 

Here we present an outline of the Lienard's Theorem for the sake of completeness. 
This helps us introducing necessary notations which will be used subsequently. 

Theorem 1 The equation has a unique periodic solution if 
(i) f and g are continuous; 

(ii) F and g (x) are odd functions with g (x) > for x > 0; 
(Hi) F is zero only at x = 0, x = a, x = —a for some a > 0; 
(iv) F (x) — > oo as x — > oo monotonically for x > a. 
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A Brief Sketch of the Proof. The general shape of the path can be obtained 
from the following observations. 

(a) Because of the symmetry of the systemQ under (x, y) — > (— ar, — y) any 
periodic orbit is symmetric about the origin. 



(6) The slope of a phase path is given by 



dy 



-9 0) 



Thus, a phase path is horizontal if 



dx y — F (x) 

■ r d y 



(3) 



dx 



0, i.e. if g (x) = 0, i.e. if x = 



(by (ii) above). Similarly, a phase path is vertical on the curve y = F(x). 
Above the curve y — F (x) we have x > and below x < 0. Moreover, y < 
for x > and y > for x < 0. 



LIMIT CYCLE y 



y=F(x) 




Figure 1: Orbits of the Lienard System 

A path YY'Y" (Figure [I]) is closed iff Y and Y" coincide, which means by 
symmetry (a) 

OY = OY' . (4) 



This is equivalent to 

Vyqy' = 

where for a typical path YQY' in Figure [2] 



Vyqy' = Vy' — Vy = J dv = J Fdy. 

YQY' YQY' 



(5) 



(6) 



and 



v(x,y) = 



g (u) du + -y 



(7) 




Figure 2: Typical paths for the Lienard Theorem 

Writing 

Vyqy' — Vyb + Vbqb' + Vb'y' 

where BB' is a line parallel to the y— axis and passing through the point (0, a) 
when the function F changes its sign from negative to positive, one then proves 
that 

(A) As Q moves out of the point A (0, a) along the curve AC, the potentials 
Vyb + Vb'Y' is positive and monotone decreasing. 

(B) As Q moves out of the point A (0, a) along the curve AC, Vbqb' is monotone 
decreasing. 

(C) From (A) and (B) it follows that Vyqy 1 is monotone decreasing to the right 
of the point A, (Figure [2l. 
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(D) The quantity Vbqb 1 tends to — oo as the paths moves away to infinity. 

(E) From (C) and (£>), it follows that the quantity Vyqy' is monotone decreasing 
to — oo, at the right of the point A (Figure [5]). 

(F) Vyqy' > when the point Q is at A or to the left of the point A. 

It thus follows from (E) and (F) that Vyqy 1 is monotone decreasing continuous 
function which changes its sign from positive to negative as the point Q moves 
out of A (a, 0) along the curve. As a result, Vyqy 1 will vanish once and only once. 
Thus, there is one and only one closed path and the proof is complete. □ 

Remark 1 The unique limit cycle in the above theorem is simple in the sense 
that no (differentiable) perturbation satisfying the conditions (i)-(iv) can bifurcate 
the limit cycle into two or more number of limit cycles. 

Remark 2 The condition (E) enables us to conclude that once Vyqy 1 becomes 
negative, it can never be positive as Q moves to infinity through the curve of 
F (x). This observation helps us to deduce the existence of a unique limit cycle. 
However, we see that the existence of the limit cycle is indeed ensured only if 
Vyqy 1 becomes negative from positive i.e., if there is a change in sign of Vyqy 1 - 
Further, the unique value of x for which Vyqy' = gives the amplitude of the 
limit cycle. Accordingly, if Vyqy' becomes negative as Q moves out from origin 
through the curve of F (x) then we get a limit cycle. This observation actually gives 
one with a possibility of weakening the conditions of the classical theorem, so as to 
accommodate a larger class of functions F (x) but still having a unique limit cycle. 
Theorem [2] is one such realizations of a stronger version of the classical theorem, 
which shows that the existence of the (unique) limit cycle actually depends on 
the local monotonicity of F (x) on a bounded interval containing the point where 
Vyqy 1 vanishes. A limit cycle can indeed be realized even when F (x) is bounded 
as \x\ — » oo (c.f. Example [l]) 

If it happens further that Vyqy 1 becomes positive from negative once more, then 
also, by an analogous argument as above we can get a point Q on the curve F (x), 
through which another limit cycle must pass. To prove this result we consider a 
function F (x) (in section |4| which is monotonically increasing to the right of the 
point A for a sufficiently large value of x and then it becomes decreasing for some 
subsequent values of x and ultimately become negative. The proof depends on an 
efficient estimate of the amplitude of the first limit cycle. 

3 Extension of the Classical Theorem and an Es- 
timate of the Amplitude 

Let, (a, F (a)) be the coordinate of Q, as shown in Figure [2] and let a = a be the 
amplitude of the limit cycle of Theorem [T] It is well known that determining the 
exact value of limit cycle of the Lienard system is a relatively difficult problem 
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[[TT], [H]]. We now find an estimate of a, for which the corresponding Vyqy' just 
become negative from positive. Since, Vyqy' is a monotone decreasing continuous 
function as the point moves out of the point A (a, 0) along the curve, without any 
loss of generality we can say Vyqy 1 can just become negative from positive if at 
least one of the following two cases hold, viz., 

(*) V YQ = but V QY > < 

(ii) V QY ' = but V YQ < o. 
The third possibility Vyq < and Vqy 1 < can be reduced to either of the above 
two cases by monotonicity and continuity of Vyqy 1 , i-e. by taking an a closer to a, 
(i.e., a — > a + ) either one of Vyq and Vqy> can be made to vanish. Similarly, the 
possibility that either one of Vyq and Vqy> is positive while their sum is negative, 
can also be eliminated by choosing a far from a (a > a). 

Case (i) 

Here Vyq = is possible if 

V YB + V bq = (8) 

In step (A) of the proof of Theorem [l] it has been proved that Vyb > 0. 
We are now going to show that Vbq < 0. 

dy 

On the path BQ, we have F (x) > and — = y = —g (x) < 0. 
Therefore, 

V BQ = J Fd y= f F %* 

BQ BQ 

= - J F{x{t))g(x{t))dt<Q. (9) 

BQ 

Thus, we can say that © is true if \Vyb\ = \Vbq\ 

i.e., if \v(a,y+ (a)) - v (0,y+ (0))| = \v (a, F (a)) — v (a, y + (o))| 



where y + (0) = OY (Figure |2]). It is possible if 

l a 

C g («) du + \y 2 + (a) - \y\ (0) - - fg (u) du - \f 2 (a) + \y\ (a) 



o 

So, we have 



G{a) = \yl{Q)-\F*{a) (10) 



where G (x) — g (u) du 
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Let a = a' be a root of (flU]) (existence of which is assured by construction) so 
that 

G{a') = \yl{Q)-\F*{a'). (11) 



Case (m) 

Here, Vqy' = is possible if 



V QB > + V B 'y> = (12) 



In step (A) of the proof of the Theorem [I] it is proved that Vb'Y' > 0. 

Proceeding analogous to case (i) one establishes that Vqb' < and conse- 



quently ( 12 ) is true provided 



G(a) = ^ (O)-V(a) (13) 
where OY' = -y_ (0) , y_ (0) < (Figure^. If a = a" be a root of (Q3J} we have 

G{a")= l -yl{Q)- l -FHa"). (14) 

It now follows that if we take 

a = max{a',a"} (15) 

then for any value of a > a, Vyqv < since, Vyqy' = Vyq + Vqy>- Thus the 
function F should be monotonic increasing in the interval a < x < a. Notice 
that the classical Lienard theorem already ensures the existence of such an a. In 
the light of the above discussion we can now extend the classical Lienard theorem 
by weakening the unbounded nature of the function F as stated in the following 
theorem and cover a more large class of functions. 

Theorem 2 The equation has a unique limit cycle if 

(i) f and g are continuous in (—d, d) for sufficiently large d; 
(ii) F and g are odd functions with g (x) > for x > 0; 
(Hi) F is zero only at x — ; x = a, x = —a for some a where < a < d; 



(iv) 3 a number a defined by (15) such that F is monotonic increasing in 



a < x < a and nondecreasing in a < x < d. 

The existence of a ensures a sign change in Vyqy' whereby we get the existence 
of a unique limit cycle in the finite phase plane. The remaining part of the proof 
of this theorem remain same as that of classical Lienard theorem. Also, from 
the above discussion and the proof of classical Lienard theorem it follows that 
Vyqy' does not change its sign any more if the function F is simply monotone 
nondecreasing in a < x < oo. In such cases the function F can even be bounded 
and even attain a constant value as x — > oo, but still we get a unique limit cycle 



8 



for such a bounded Lienard system. Thus, we can indeed cover a larger class of 
functions than those covered by the Lienard theorem (c.f. [l}. 

In the beginning of the proof of Theorem [I] we observed that above the curve 
y = F (x) we have x > and below x < 0. So, the ir-coordinate of a point on 
a limit cycle will achieve its maximum absolute value on the curve y = F(x). 
Therefore, the amplitude of a limit cycle for the Lienard system is the abscissa 
of the point Q lying on the curve y = F(x). Since, for a limit cycle we have 
Vyqy' = 0, so by the construction of a it follows that it is an efficient upper 
estimate of the amplitude of the limit cycle. In the following example we find the 
values of a for the well known Van der Pol equation against different values of \x 
and compare them with the results obtained in [TT] and [T3]. This also gives an 
example of a bounded Van der Pol equation having same amplitude as that of the 
standard Van der Pol equation. 

Example 1 Here, in the following table we present estimates of the amplitude 

fx 3 \ 

of the limit cycle for Van der Pol equation in which F (x) — /i I — — x J and 

g (x) = x for different values of fi. It is clear that our estimates are reasonably 
close to the exact (numerically computed) values (as reported in jH]). Our values 
also appear to be much better than the upper bound 2.3233 of [11] (the estimated 
values of [14] are valid only for small fi). From our numerical estimates it follows 
that although the estimated values of the amplitude seem to vary irregularly for 
the moderately large values of [i £ [0, 10], these are nevertheless bounded above 
by 2.05. 



t 1 


y+ (0) and y^ (0) 


a 


0.1 


2.00117 


2.0000586437166383 


0.2 


2.007076 


2.002540101136999 


0.3 


2.015912 


2.0054678254782505 


0.4 


2.028253 


2.0091503375996034 


0.5 


2.044065 


2.013278539452526 


1 


2.1727135 


2.0327736318429275 


1.5 


2.3710897 


2.0436704679281523 


2 


2.6149725 


2.04739132291152 


2.5 


2.8844602 


2.047213463900291 


3 


3.1687156 


2.045311842105752 


3.5 


3.462322 


2.0427848260891426 


4.5 


4.06701715 


2.037557405718347 


5 


4.3752293 


2.035154629371522 


10 


7.5528123 


2.020095969119061 
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We get the same result if we consider the bounded function 
/' ' \ ( 

4,76/1 
sin (0.6) 



F(x)= < 



0T 
3 

3 

(2.4) 3 
3 

(2.4) 3 



- 2.4 
2.4 I 



sin (0.6) 



a; G (-2.4,2.4) 
(cos (0.6) - cos (a: - 3)) a; G (-3, -2.4] U [2.4, 3) 
(cos (0.6) -1) x G (—oo, —3] U [3, oo) 



in the whole phase plane or the following function in finite phase plane. 



F(x) 



(2.4)< 



-2.4 - 



4.76// 
sin (0.6) 



x G (-2.4,2.4) 
(cos (0.6) - cos (x - 3)) x G (-3, -2.4] U [2.4, 3) 



This also tells that the value of d need not by too large. A limit cycle is assured 
even for a moderately large d. 

It follows from this example that the amplitude of the unique limit cycle is 
independent of the asymptotic behaviour of F (x) as x — > oo. Indeed, the ampli- 
tude corresponds to the point Q (a, F (a)) on the limit cycle for which the integral 
Vyqy' — J F dy vanishes and clearly depends on the form of F (x) in the fi- 

YQY> 

nite interval (0, a). To the best of author's knowledge, this result apparently is 
not recorded clearly in the literature. We therefore state this observation as the 
following corollary. 

Corollary 2.1 Amplitude of the unique limit cycle of the Lienard system |2]) is 
independent of the asymptotic behaviour of F (x) as \x\ — > oo. 



4 The New Theorem 

By the observations discussed in last section it is clear that, in the interval a < 
a < a we obtain the limit cycle of Theorems [l] and [2] Moreover, because of the 
condition (iv) we see that the limit cycle remains unique. However, if the function 
does not satisfy this condition, then the limit cycle may not be unique. We now 
present our new theorem. 

Theorem 3 Let f and g be two functions satisfying the following properties, 
(i) f and g are continuous; 

(ii) F and g are odd functions and g (x) > for x > 0.; 

(Hi) F has positive simple zeros only at x — a\, x = 02 for some a\ > and 



some 02 > a, a being defined by (15) and a < L, where L is the first 
local maxima of F (x) in [a± , 02] ; 
(iv) F is monotonic increasing in ai < x < a and F (x) — > —00 as x — > 00 
monotonically for x > 02] 
Then the equation ([T]) has exactly two limit cycles around the origin. 
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Proof: We can get exactly the same observations as we get in observations 
(a) and (b) in beginning of the proof of Theorem IT] 




Figure 3: 

By the observations in section [3] we can ensure the existence of inner limit 
cycle. So, we shall now prove the existence of one more limit cycle by showing 
that 

OY = OY' 

once more when x > a. To prove the result we shall consider the function v {x, y) 
as in ([7]), and write, 

V Y qy> = V YX + V XB + Vbqb- + Vb-x> + Vx>y> (16) 

where, XX' is a line parallel to the y— axis passing through the point (0, a±) where 
the function F changes its sign from negative to positive and BB' is a line parallel 
to the y— axis passing through the point (0,02) where the function F changes its 
sign from positive to negative. The proof is carried out through the steps (A) to 
(F) below. Here we refer to the Figure [3j 
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Step (A) : As Q moves out from A2 along A2C, Vyx + Vx'Y' is positive 
and monotonic decreasing. 

We choose two points Q (a, F (a)) and Qi (ax,F («i)) on the curve of F (x) where 
«i > a. Let YQY' and YiQ\Y{ be two paths through Q and Qi respectively. On 
the segments YX and Y\X\ we have 

y > 0, F (x) < and y - F (x) > 0. 

Now, 

(y-F(x)) YX < (y-F(x))y iXi 
> 



^y-F(x)J YX \y-F(x) /YiXi 
Since g (x) > for x > 0, we have 



-9(x) \ ( -g(x) 



y-F(x)J YX \y-F(x) /YiX 



So by (§ we get 



Therefore, 



^ <m <» ,i7) 



Using (17) we get 



J l^r j£ \^ diX J -^r j£ 



V YX = J Fdy= J (-F) )'!■'■ 

YX YX 



Y1X1 

Since, F and dy are positive on Y\X\ we have 



VVx > J Fdy = V YlXl > 0. (18) 



YlXi 

Next, on the segments X'Y' and X[Y{ we have 

y < 0, F (at) < and y - F (x) < 0. 



Now, 



So, by (§ 



-9(x) \ ( -g(x) 



y-F(x)J x , Y , \y ~ F (x) J xlY , 



>(*'] >0. (19) 



dx J j 1 ^ 1 ^^ 1 V dx 
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Therefore, 



V X nr>= J Fdy= J {-F)^dx 



Using ( 19 1 we get 



X'Y' Y'X' 



V X , Y ,> J (-F)f x d X 



Since, F and dy are negative on X[Y{ we have 

Vx'Y- > J Fdy = V X [Yl > 0. 



(20) 



From (18) and (20 1 we have 

V Y X + V X 'Y' > V YlXl + V x[ yi > 0. 

Therefore, Vyx+Vx'Y' is positive and monotonic decreasing as the point Q moves 
out from A-i along A-iC. 

Step (B) : As Q moves out from A 2 along A 2 C, V X b + Vb'x 1 ls negative 
and monotonic increasing. 

On the segments XB and X%Bi we have 

y > 0, F (x) < and y - F (x) > 0. 

Now, 

(y-F(x)) XB < ( y -F(x)) XiBi 



~9 (%) 



< 



y-F(x)J XB \y-F(x)J XiBi 



-g{x) 



So, by ([3]) we get 



< 



XB 



< 



XtBt 



Therefore, 



V; 



XB 



Fdy = .' I '^d.r 



XB XB 



Using (21 ) we get 



V XB < J F^dx = J Fdy = V XlBl < 



since, F > and dy < on XB and XxB±, 



(21) 



(22) 
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Next, on the segments B'X' and B' 1 X' l we have 

y < 0, F (x) > and y - F (x) < 0. 



Now, 



(y-F(x)) B , x , > (y-F(x)) 



B'X' 



'9 (x) 



y-F{x) 



> 



B'X' 



-9 (x) 
y-F (x) 



Using (|]) we get 



> 



B'X' 



> 0. 



B[X[ 



Therefore, 



V, 



B'X' — 



F,hj=- I F^dx= ( F(-^\d, 



B'X' 



X'B' 



X'B' 



So, by ( 23 ) we have 

V B 'x' < 



Fdy = V B , x , <0. 



(23) 



(24) 



since, F > and dy < on B'X' and B[X[. 
From ((22| and pi| we have 



VxB + VB'X' < VjfiBi + Vb^xj < 0. 

Therefore, Vxb+Vb'X' is negative and monotonic increasing as the point Q moves 
out from A 2 along A 2 C. 

Step (C) : As Q moves out from A 2 along A 2 C, Vbqb' is positive and 
monotonic increasing and tends to +oo as the path recedes to infinity. 

On BQB' and BiQiB[, we have F (x) < 0. We draw BH X and B'H[ parallel to 
x— axis. 



Therefore, 



Vb 1 q 1 b' 1 = j Fdy 
b 1 q 1 b[ 



(-F) dy 



B' l Q 1 B 1 



> J (-F)dy 
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since, F (x) < and dy > for points on H[QiHi. Again since, 



F{x)] BIQB > F(x)] H , QiHi 

for same value of y we get 



V BlQl Bi> J (~F)dy> J (-F)dy 

i^Qi-Hi B'QB 

Fdy 



BQB' 



=> V BlQlB { > Vbqb-- (25) 

Next, let S be a point on the curve of F (x), to the right of A 2 , and let BQB' 
be an arbitrary path, with Q to the right of S. The straight line PNSP' is parallel 
to the y— axis. Then, 



V B qb>= J F(x)dy 



BQB' 

{-F(x))dy 

B'QB 

> J (-F(x))dy (26) 
P'QP 

since, (-F (x)) > and dy > along B'QB. Now by condition (iv) of this theorem 
it follows that F is monotonic decreasing for x > a 2 and so we have \F (x)\ > NS 
on P'QP and since further F (x) < on P'QP so this implies — F (x) > NS on 
P'QP. Again, PP' > NP'. Thus we get 

V B qb' > J NSdy = NS J dy = NS- PP' > NS ■ NP' . 

P'QP P'QP 

But as Q goes to infinity towards the right , NP' — > oo. Hence, we can conclude 
that Vbqb 1 is positive and monotonic increasing and tends to +oo as the paths 
recede to infinity. 

Step (D) : 

From steps (A) and (B) it follows that the quantities Vyx+Vx'Y' and Vxb+Vb'X' 



are bounded quantities. Thus by ( 16 ) and by step (C) it follows that Vyqv is 



monotonic increasing to +oo to the right of A 2 . 
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Step (E) 



By the construction of a it is clear that Vyqv < in a < x < a 2 i.e., to the left of 
A 2 . Again from step (D) we conclude that Vyqv ultimately becomes positive as 
Q moves out of A 2 along the curve of F (x). Therefore, by the same reason given 
in conclusion of the Theorem [TJ it follows that there is one and only one path in 
the region x > a such that 

Vyqv = 0. 

Also, by ([5]) and the symmetry of the path it is clear that the path is closed. 
Step (F) : 

By the construction of a and by step (E) it is clear that equation (JIJ has exactly 
two limit cycles around the origin, the second limit cycle surrounds the first one. 
This completes the proof of the Theorem [3] 

Remark 3 It also follows from the proof that both the limit cycles are simple (c.f., 
Remark [TJ that neither can bifurcate under any small C 1 perturbation satisfying 
the conditions of the theorem 

Remark 4 One cannot assume that Vyqy' < if a > L. We give a counter 
example below. 

Remark 5 It is well known that two consecutive limit cycle cannot both be stable 
(unstable) . Because of our choice of the function F (x) (negative and monotone 
decreasing at the right of and near the origin and infinity) , the inner limit cycle is 
stable and outer limit cycle is unstable (in reverse to those of reference [jS], [IS]])- 

The existence of exactly N limit cycles is established by an easy extension of 
the above proof [19]. We state the theorem as follows. A brief outline of its proof 
is given in the Appendix. 

Theorem 4 Let f and g be two functions satisfying the following properties, 
(i) f and g are continuous; 

(ii) F and g are odd functions and g (x) > for x > 0.; 

(Hi) F has N number of positive simple zeros only at x — o» } i = 1,2, . . . ,N 
where < a\ < a 2 < . . . < ajy such that in each interval Ii — [at, dj+i], 
i = 1, 2, . . . , N — I, there exists cti, satisfying properties given by \15\ , 
such that c\i < Li where Li is the unique extremum in Ii, 
i = 1, .. .,N — 2 and Ln—\, the first local extremum in [ajv-i,ajv]- 
(iv) F is monotonic in ai < x < a t V i and \F (x)| — » 00 as x — > 00 
monotonically for x > . 
Then the equation (T7| has exactly N limit cycles around the origin, all are simple. 

Remark 6 The conditions (1) and (iv) of Theorem [3] and Theorem [4] may be 
weakened following Theorem [2] For instance, the condition (iv) of Theorem [3] 
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may be restated as 

(iv) F is monotonic increasing in a\ < x < a = a\ (say) and 3 a number oti 



given by ( 15 ) such that F is monotonic decreasing in 02 < x < a.2 and 



nonincreasing in oli < x < d. 

5 Examples 

It is shown in section 3.3 of [TS] that the limit cycles of the autonomous system 



x = y 
-x — fih (x, x) 



(27) 



are asymptotic to the circle x 2 + y 2 = r 2 as fi — > where the values of r are the 
roots of the equation 

/-2tt 

$ (r) := / h (r sinu, r cosu) cosu du = 0. (28) 
^0 

We note that this is not the Lienard system. It is the canonical phase plane for 
Lienard equation. The phase diagram of above system and the Lienard system, 
however, should be similar. Here we take fj, = 0.1, h (x, x) = (—4 + 75a; 2 — 50/cx 4 ) x, 
k 7^ so that 



f(x)=fi (-4 + 75.x 2 - 50fca; 4 ) = -0.4 + 7.5x 2 - 5k 
F (x) = -OAx + 2.5a; 3 - kx 5 



and 

f-271 



^ ( r ) = I ( — 4 + 75r 2 sin 2 u — 50fcr 4 sin 4 u) r cos u ■ cos u du 

3 

1 



— 7rr (25fcr 4 - 75r 2 + 16) 



Therefore, ( 28 1 reduces to 

1 



-~wr [25kr 4 - 75r 2 + 16) = 



giving 



r 2 

10k 



— (l5 ± V225 - 64fc) 

UK \ / 



So, (28 1 has real and distinct roots if 225 > 64fc i.e., if k < 3.515625, real and 
repeated if k = 3.515625 and imaginary if k > 3.515625. Therefore it follows 
that we will get two distinct limit cycles, which are asymptotic to the circles 
corresponding to the above two distinct values of r if k < 3.515625. Similarly, 
we will get only one limit cycle when k = 3.515625 and no limit cycle when k > 
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3.515625. It can be verified that the system undergoes a saddle node bifurcation 
at k — 3.515625. 

We note that the point (a, F (a)) on the limit cycle in the Lienard plane gets 
transformed to the point (—a, 0) lying on the almost circular limit cycle of radius 
r > r! (in the canonical phase plane) where r\ = ^ (15 — \/225 — 64fc) (rj 
corresponds to the first limit cycle) under the transformation 

x = —u, y = —v + F (u) 

(u,v) and (x, y) being the corresponding points in Lienard plane and canonical 
phase plane with / an even function. We thus have 



> \l^r f 15 — V225 - 64fc 



lOfc 

We now present the phase diagrams of the above systems in Lienard plane in 
the following examples for different values of k. These examples justify our new 
theorem. We use Mathematica 5.1 in constructing the examples. 

Example 2 Here we consider the autonomous system ^ with k = 3.65, / (x) — 
-0.4 + 7.5a; 2 - 5fcr 4 g (x) = x and F (x) = -0.4a: + 2.5a; 3 - kx 5 . The phase 



diagram in Lienard plane is shown in Figure 4(a) which does not have any limit 
cycle. Again we take k = 3.57 in the above system. The corresponding phase 
diagram is shown in Figure |4(fe)| This phase diagram also does not contain any 
limit cycle, but we see that the path is concentrating in a certain circular region. 

Example 3 Here we consider the autonomous system |2]) discussed above with 
k = 3.515625 so that a x ~ 0.49307, a 2 ~ 0.68410, y + (0) and y_ (0) both are 
approximately equal to 0.652287 and a ~ 0.65204. Let L% be the point of minima 
of F in (0, ai) and Li be the point of maxima of F in (01,02). Here F (x) is 
increasing in 

(-i 2 ,ii)U(ii,i 2 ) 

and decreasing in 

(-oo,-L 2 )U(-L 1 , J Li)U(L 2 ,oo) 
where, Li ~ 0.24997 and L 2 ~ 0.60348. In this case we obtain only one limit 



cycle as shown in Figure 4(c) Here, F is not monotone increasing throughout the 
interval a\ < x < a, violating the condition (iv) of Theorem [3] (since a > L 2 ). 

Example 4 We now take k = 3.5. Here, y + (0) and y_ (0) are approximately 



equal to 0.624499. The equations (10 1 and (131 both reduce to 



Y = \ (0-624499) 2 - i (-0.4a; + 2.5x 3 - 3.5a; 5 ) 2 

having real roots x = ±0.62393 so that a ~ 0.62393. Here, a\ = 0.4919 and 
a 2 = 0.68725 showing that ai < a < a 2 . Next, L Y ~ 0.24985, L 2 ~ 0.60510. Here 
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(a) 



(b) 




(c) 



(d) 



Figure 4: The phase diagram of the system ^ in Lienard plane 

(a) with k = 3.65, and center as a repelling node in Example [2j 

(b) with k = 3.57 in Example [5J 

(c) with k — 3.515625, and one limit cycle in Example |3j 

(d) with k — 3.5, and two limit cycles in Example [4j 
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all the conditions of Theorem [3] are satisfied except condition (iv). However we still 
get two limit cycles as shown in Figure [4 ( d) | drawn in Lienard plane. This example 
and the above example show that the conditions of Theorem [3] are sufficient but 
not necessary. 

Example 5 Finally we take k = 3. Here, y+ (0) and y_ (0) are approximately 
equal to 0.5552 and a ~ 0.55324. Here, a x = 0.46473 and a 2 = 0.78572 showing 
that a\ < a < a%. Next, L\ ~ 0.24638, L 2 — 0.66279. Here, all the conditions of 
Theorem [3] are satisfied and so we get exactly two limit cycles. The phase diagram 
in Lienard plane is shown in Figure [5j 




Figure 5: The phase diagram of the system @ in Lienard plane with k — 3, and 
two limit cycles in Example [5] 



Remark 7 Although in the above examples the value of /j, is sufficiently small 
(so as to satisfy the amplitude analysis of [TS]) our theorem should be applicable 
for large values of More detailed bifurcation analysis in the (/z, k) parametric 
plane will be considered separately. 



Example 6 We now consider the function 
-0.1 sin (IOto) 



F-i 



0.01 



0.15 



0.01 

0.02099503719021- 
2^/0.1 (x - 



< x < 0.15 
0.15 < x < 0.15 

x > 0.15 



T01 



1 



x 101 



0.2395037190209989) 
Then we have oi = 0.1, a 2 = 0.25052350868645645, L = 0.15. Here, 

/ (£) = < / (0.2395037190209989) = 0.2006848039831627 > / (a 2 ) = 0.9526060763219791, 
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though 

L < 0.2395037190209989 < a 2 

showing that the function is not monotone nonincreasing in [L, a 2 ] and so it 
does not satisfy the condition (3) of Theorem 5.1 in chapter 4 in the book 0]. 
However for the inner limit cycle we have y+ (0) = y_ (0) = 0.12238318. So, 
a = 0.12221435874426823 < L satisfying the conditions of Theorem [3| This ex- 
ample clearly shows that the Theorem [3] covers a larger class of functions than 
those covered by Theorem 5.1 in chapter 4 in [4]. The function F alongwith two 
limit cycles are shown in Figure [6j 




Figure 6: The phase diagram of (J5| in Lienard plane for Example [6j 



Example 7 We now consider a different problem. Here we define 

- 0.048989794 \ 2 
O05 J 



where 



0.005 - 0.025W 1 - 



F+ (x) = { 



-0.0008137888130718 + 0.01 W 1 - 



0.14781375 V 
0.05 J 



< x < ai 
a% < x < ci2 



x _ n 29746094 \ " 

0.()0(l91(i.sll6(l(i400276.-) - 0.015, /I - | — J «., - . ,,■ .. „■, 



-0.0003265987749816556 + 0.04Va; - 0.3972073012751128 x > a 3 

ax = 0.097979588 
a 2 = 0.197647912 
and a 3 = 0.397273968. 
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and 

F+ (x) x>0 



F(x) = 



-F + (-x) x < 



The function F+ (x) is obtained by matching three ellipses and a parabola succes- 
sively in the intervals (0, a\) , (cii, a 2 ) , {a 2 , 03) , and (03, 00) such that 

F+ (a, + 0)=F + (a, - 0) and F' + (a, + 0) = F' + (a, - 0) , (29) 

where etj's are zeros of F + . The unique extremum of F in (0, ai), (ai, a 2 ), (02, 013) 
are respectively 

L a = 0.048989794 
Li = 0.14781375 
and L 2 = 0.29746094. 

We obtain three limit cycles which meet the positive y-axis at the points (0, y\ (0)), 
(0,j/ 2 (0)), (0,y 3 (0)) where 

2/i (0) = 0.1332869 
2/2 (0) = 0.212146685 
and 2/3 (0) = 0.4630114. 



The matching conditions (29) are used to make F € C 1 (i?) with accuracy level 
O (l0~ 7 ). This function is constructed in a trial and error method and numer- 
ical data with large significant digits arise in this fashion. Examples with lower 
significant digits and lower and higher accuracy are possible in principle. Here 
we get at = 0.133002186 and a 2 = 0.21203506657. The function F satisfies all 
the conditions of Theorem [4] (for example on < Li etc.) and so the existence of 
the above three limit cycles are ensured by this theorem, the proof of which is 
presented separately [19] . However, the function F is defined in such a manner 
that |F(Lo)| > |F (L2)| implying that j3 2 mentioned in Theorem 1 of [17! or in 
Theorem 7.12, chapter 4 of the book [4] does not exist and hence these theorems 
are not applicable for the corresponding Lienard system. The limit cycles of the 
Lienard system in Lienard plane and the graph of the function F have been shown 
separately in the Figure [7j To conclude, Theorem 1 of [17 or Theorem 7.12, 
Chapter 4 in the book [3] fail to predict the existence of the exact number of limit 
cycles for the above function F (x). 



6 Concluding Remarks 

Many interesting new results have been proved on the existence of an exact number 
of multiple limit cycles [[8]. [15] ■ [IT]] in the recent past. Odani has proved a 
sufficient condition in 8J using a choice function cf>k which can be exploited to 
obtain better estimates of amplitudes of the limit cycles. We used a straight 
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(a) (b) 



Figure 7: 

(a) The phase diagram of the system ^ in Lienard plane with three limit cycles. 
(6) Graph of the function F in Example [7j 



forward method depending on the geometry of phase diagram. We have proved 
a similar result with more general class of functions F (x) by a simpler method. 
In the present approach a strict monotonicity of F (x) is required only in the 
intervals ai < x < a and x > 02- Consequently, F (x) can accommodate "small 
scale" oscillations in the interval a < x < 0,2- Odani, for instance, considered 
an F (x) which is not only C 1 but also has a unique extremum in the interval 
<Xi < x < 0,2- Further, the theorem is valid for a more general class of the function 
g (x). Odani's theorem however is valid only for g (x) = x. An interesting problem 
will be to establish the relation between a of our approach and the function </>£.. 
We note that on corresponds to the amplitude of the limit cycles. Our estimates of 
amplitude of the limit cycle of the Van der pol equation constitute an improvement 
over those available in the literature [[H], [II]]. Examples [6] & [7j on the other 
hand, show the difference between the present theorem and those of [4] and [T7| . 
The calculations are accurate upto the accuracy level O (10~ 7 ). The existence of 
limit cycles in a Lienard system allowing discontinuity (see, for instance, |18j ) is 
an interesting problem for further study. Determining the shape of the limit cycles 
is also left for future investigations. 

Before closing we note that the value a, in general, is a function of the pa- 
rameters of F (x) in the parametric space. For instance, in Examples [2] — 1~5| a 
is a function of the parameters fi and k. The study of the variation of a in the 
parametric space seems to offer interesting insights into the bifurcation and re- 
lated issues of the multiple limit cycles in a Lienard system. The relationship 



23 



with Poincare's return map also needs to be studied. We wish to investigate these 
problems in future. 

7 Appendix 

Here we present a brief outline of the proof of Theorem [4] which is given in detail 
in |19j . The theorem is proved by an induction method which is dependent on the 
non-trivial initial hypotheses that the the result holds for N = 1 and N = 2. 




Figure 8: 

We shall prove the theorem by showing the result that each limit cycle intersects 
the a;— axis at a point lying in the open interval (cSj, Qj+i], i = 0, 1, 2, . . . , N — 1, 
where <5o = Lo is the local minima of F (x) in [0,eti]. By Lienard theorem and 
Theorem [3] it follows that the result is true for TV = 1 and N = 2. We shall now 
prove the theorem by induction. We assume that the theorem is true for N = n— 1 
and we shall prove that it is true for N ~ n. We prove the theorem by taking n 
as an odd +ve integer so that (n — 1) is even. The case for which n is even can 
similarly be proved and so is omitted. It can be shown that pQ, Vyqy' changes 
its sign from +ve to —we as Q moves out of A\ (a l5 0) along the curve y = F (x) 
and hence vanishes there due to its continuity and generates the first limit cycle 
around the origin. Next, in Theorem [3] we see Vyqy 1 again changes its sign from 
—ve to +ve and generates the second limit cycle around the first. Also, we see 
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that for existence of second limit cycle we need the existence of the point a, which 
we denote here as a.\. 

Since by induction hypothesis the theorem is true for N = n — 1, so it follows 
that in each and every interval {a k , <5fc+i], k — 0, 1, 2, . . . , n — 2 the system ^ 
has a limit cycle and the outermost limit cycle cuts the x— axis somewhere in 
(a„_i,oo). Also Vyqv changes its sign alternately as the point Q moves out of 
di's, i = 1, 2, . . . , n — 1. Since (n — 1) is even, it follows that Vyqy< changes its sign 
from +ve to —ve as Q moves out of a„_ 2 along the curve y = F (x). Since there is 
only one limit cycle in the region (a n -i, oo), so it is clear that Vyqy' must change 
its sign from —ve to +ve once and only once as Q moves out of A n _\ (a n _ 1; 0) 
along the curve y = F (x). Also it follows that once Vyqv becomes +ve, it can 
not vanish further, otherwise we would get one more limit cycle, contradicting the 
hypothesis so that total number of limit cycle becomes n. We now try to find an 
estimate of a for which Vyqv vanishes for the last time. 

We shall now prove that the result is true for N — n and so we assume that 
all the hypotheses or conditions of this theorem are true for N = n. So, we get 
one more point a n and another root a n , ensuring the fact that Vyqy' vanishes as 
Q moves out of A„_i through the curve y = F (x), thus accommodating a unique 
limit cycle in the interval (a n -i,a n \. 

By the result discussed so far it follows that Vyqy' > when a lies in certain 
suitable small right neighbourhood of a n _i. We shall prove that Vyqy' ultimately 
becomes —ve and remains —ve as Q moves out of A n (a n ,0) along the curve 
y = F (x) generating the unique limit cycle and hence proving the required result 
for N = n. 

We draw straight line segments XkX' k , k = 1, 2, 3, . . . , n, passing through Al- 
and parallel to y-axis as shown in Figure [8] For convenience, we shall call the 
points X n , X' n , Y, Y' as B,B' , Xq, X' respectively. We write the curves 
Y k = X k - 1 X k , r k = X' k X' k _ v k = 1,2,3,..., n 

so that 

n n n 

yqy' = x qx' = E r fe + x n QX' n + £ r' fc = E (T fc + ri) + bqb< 

k=l ft=l k=l 

and 

Vyqv = E (Vv„ + V r > k ) + V BQB , . (30) 



which is used in place of the function in (16). The rest of the proof are analogous 
to that of Theorem [3] and proved separately in [19] . 
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